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ABSTRACT 


A series of average transmission and average self-indication 
ratio measurements have been performed in order to investigate 
the temperature dependence of the resonance self- shielding effect 
in the unresolved resonance region of depleted uranium and tanta- 
lum. The measurements were carried out at 77°K ) 295°K and 
e»1000°K with sample thicknesses varying from ~0 . 1 to ~1.0 mean 
free path. The average resonance parameters as well as the 
temperature dependence have been determined by' using an analyti- 
cal model which directly integrates over the 'resonance parameter 
distribution functions. 
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INTRODUCTION 

The main contribution of the Doppler effect in a fast breeder 
reactor comes from "the unresolved resonance region"^where the 
detailed resonance structure has not been, and possibly cannot be, 
experimentally determined. The temperature variation of the neu- 
tron f^ux inside the reactor must be predicted for safe reactor 
operation The prediction of this temperature -dependent varia- 
tion of neutron flux follows from: 

1) the determination of the average resonance parameters 
with which the cross section and resonance structure can be repro- 
duced (or at least simulated), 

2) sufficient knowledge about the temperature dependence of 
the resonance self- shielding effect based on the above average 
parameters. 

In order to predict the cross section in the unresolved 
region it is necessary to assume that the resonance parameters 
are governed by some specified statistical distribution functions 
(such as the Wigner level spacing distribution^ and the Porter- 
Thomas neutron width distribution 3 ). Furthermore it is necessary 
to assume that the average resonance parameters determined in the 
resolved region can be extrapolated into the unresolved region. 

One of the main purposes of this investigation is the extablish- 
ment of the validity in extrapolating the low energy parameters 
of Ta and depleted U into the unresolved region, using only a 
single level formalis^m. 

One of the best ways to achieve the above goal is to perform 
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a series of direct measurements on a system with good geometry 
and to obtain as many sets of data as possible by varying sample 
temperature and thickness. The average self-indication ratio, < SIR>, 
and the average transmission, <TR>, measurements were selected for 
the experimental investigation of the Doppler effect in the un- 


resolved region. Reviewing previous temperature-dependent meas- 
urements, no <SIR> measurements have been reported although the 
data from <TR> measurements^’^ are presently available. This is 
due to the difficulty arising from the unusually large background 
-to- signal ratio in the <SIR> measurements. The experiments un- 
dertaken here are the first to attempt direct measurements by 
combining both <TR> and <SIR> methods. 

These data were analyzed by using both a statistical sampling 
method^ ’ ^ and an analytical method. The results obtained by 

6 V 

using the statistical method (computer code DAISY ) have been pre- * 

g i o /OH o-v*. &L u+**£k*-* 

sented in previous publications 0 ’ ^and will not be presented 

here. Both the ideal gas model and the Lamb effective temperature 

model'*''*' were applied and compared. 


EXPERIMENTAL TECHNIQUE 

Most of the aspects of the experimental technique have been 

8 10 

described in previous reports ’ and are briefly summarized here. 
The RPI linear accelerator was used to produce neutrons for the 
time-of- flight experiment. The accelerator was operated at an 
electron energy 70 MeV, a repetition rate 480 pps and with an 
electron pulse width of 20 to 250 nanoseconds. 
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A B-Nal detector at a 28 meter flight path and a large scin- 
tillator detector at 25.69 meters were used along with a PDP-7 on- 
line computer for data acquisition. * 

The self-indication experiment along with a typical self- 
indication spectra (for the 66.1 eV resonance) are illustrated in 
Figure la and lb. These measurements were performed at 77 + ^°K 
and 295+2 K for both tantalum and depleted uranium (99.8% and 

0.2% u) and 973+5°K for depleted uranium (1073+5°K for tantalum) 
varying the sample thickness from -0.1 to -1.0 mean free path. The 
overall view of the electronics and the data acquisition system 
are shown in Figure lc. An evacuated furnace and a liquid-nitrogen 
cryostat were used to obtain high temperature samples and for low 
temperature samples, respectively. An automatic sample cycling 
system was utilized to alternate sample "into ,r and "out of" the 
neutron beam. The capture sample inside the' self-indication de- 
tector was kept at room temperature throughout- the self-indica- 
tion experiments. 

The average self-indication ratio, <SIR>, and the average 
transmission, <TR>, are defined as: 


<SIR> 


<TR> 


I,., Y Y ( E > T 2 n 2) ex P^ n i (:T t ( E >T 1 ) ]dE 

Hi 

L Y V <E ’ T2,I, 2 )dE 

Ai S AE exp[ -V't< E > T i> :,dE 


( 1 ) 


( 2 ) 


where the subscripts 1 and 2 indicate the shielding and self- 
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indication samples, respectively. The sample temperature, and 

12 

sample thickness are denoted by T and n. The capture yield, 

Y , can be approximated by for ^ <<: l- 4 

ANALYSIS 


The data were averaged over 10%-wide energy bins after dead- 
time and background corrections were made. The counting statis- 
tics (in- standard deviations) of the experimental results are 
typically •+ 1% in <TR> and + 2 % in <SIR>. Since very thin self- 
indication samples (0.00147 atom/barn for Ta and 0.00379 atom/ 
barn for depleted uranium) were used, the multiple scattering 

effects in calculating <SIR> can be neglected in the unresolved 

10 

resonance region. 

The analytical method has been developed and applied to 
interpret the data. In this method, the expressions for <TR> and 
<SIR> in Eqs . (1) and (2) are separated into two parts: the un- 
shielded term and the shielded term. Assuming that 
| a t - <a t >| <l/n^and «1, the definitions of <TR> and <SIR> in 
eqs. (1) and (2) become 


'SIR> = 


- n l < °t > 

<a > 

Y 


<a Y e 


- n i<°t • <0 t > > 


> 


(3) 


<TR> = e 


-ni < c t > <e -n l( o t - 


<a 


t». 


(4) 


where <. . .> stands for the average over the energy interval, AE. 
Doppler broadened total and radiative capture cross sections are 
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denoted by a t and a , respectively. 

The second terms in eqs. (3) and (4) are defined as the 
self-shielding factors and , and the exponential expres- 

sions appearing in r) S j R and rj TR are expanded into four terms 


^TR 




-nj^Ca - <a >) 

<e C t > 

1 + 7 n l <( ^ a t " <a t >) 2 

- \ n i < (^ t - <a t >) 3 > 


^SIR 


<o 

y 




< o t >) 


(5) 


(5) ’ 

( 6 ) 


~ <a Y > - n i <a Y ( a t - <a t >)> 

, 1 2 ^ , 2 
+ J n l <a Y ( a t " <a t >) > 

' Z n l <a y (CT t " <a t >)3> » . (6)’ 

where the linear term in eq. (5)' vanishes since <(a. - <a >)> 

t t 

- 0. 

The self-shielding factors are functions of the average 
cross sections <a y > and <a fc >, the variance term <(a fc - <o t >) 2 >, 

the triance term <(a fc - <a t >) 3 >,and the values of these respec- 
tive terms weighted by the capture cross section. Since the 
evaluation of these terms by energy integration is impossible 
without the detailed information about point-wise cross sections 
in the unresolved region, the energy integrations were replaced 
by the integrations over the corresponding statistical distribu- 
tion functions: viz. the Wigner level 
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2 

spacing distribution and the Porter-Thomas neutron width distri- 

3 

bution functions. This analytical model which directly integrates 
over the statistical distribution functions is 10 to 10 faster 
in computation time than the statistical sampling method. ^ The 
Breit-Wigner single level formula was used for s- and p-wave 
cross sections. The cross section overlapping effects of the 
same resonance sequences as well as the different resonance 
sequences were considered. The details of the analytical method 
have been presented elsewhere . For given sets of average reso- 

nance parameters , <TR> and <SIR> were calculated and compared 
with experimental values for different sample thicknesses, energy 
intervals, and sample temperatures. 

The multi-level effects which can be significant in the 
fissile nuclides and medium-weight nuclei were not considered. 
However, the experimental evidence shows that the influence of 
multi-level effects on the temperature dependence is almost 
negligible. ^ 

DEPLETED URANIUM RESULTS 
A) Average Resonance Parameters 

The average resonance parameters which best fit the experi- 
mental results (<TR> and <S1R>) in the energy range of 10 to 100 
keV are tabulated in Table I. In Figures 2a and 2b the results 
are plotted in terms of <SIR> and <TR>, and the comparison of 
these results with the ENDF/B-III evaluated parameters is given. 

In Figures 3a and 3b the average total and capture cross sec- 
tions calculated using the above parameters (denoted by grand 
average parameters in the figures) are compared with those cal- 
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culated using the local average parameters which best fit the <TR> 
and <SIR> of each 10%-wide energy bin for five different sample thick 
nesses^ p_ wave strength function, S-^, which best, fits the ex- 
perimental results ranges from 1.27x10 ^ to 2.4x10 for corres- 
ponding variations of the scattering length, R , ranging from 
9.0 to 9.6 fermis. Since the potential cross section, Op> cal- 
culated from these values of R’ (Up = 10.0 barns for R 1 = 9.0 fm 

and a =11.4 barns for 9.6 fm at 10 keV), are within the error 

P 14 15 

limitations of measured values, 5 the choice of any set of 
parameters (labeled A,B, or C in Table I) can be made only if the 
potential cross section data are improved further. The sensitiv- 
ity of <TR> and <SIR> to the changes in R* , strength functions, 
and level spacings are shown in Figures 4a, 4b and 4c. Accord- 
ing to Figure 4b, the <SIR> is more sensitive to the change in 
p-wave strength function, S-^, while the <TR> is more sensitive to 
the change in s-wave strength function, S Q . This makes the <SIR> 
experiment a powerful method for a determination of both s- and 
p-wave resonance parameters when combined with the <TR> 
experiment . 

It should be noticed that the <SIR> of the depleted uranium 
could not have been fit well by using the conventional (2J+1) 1 - 
dependence formula for the level spacing which assumes parity 
independence. The best- fit p-wave level spacing, <D^> is 
11.3 + 3.0 eV, which is approximately 60% larger than <D^> cal- 
culated from < D q > according to (2 J+l)" 1 - dependence formula for 
<D. ,>. This could indicate that there is parity dependence of 

Xj J 
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level spacing. The fit obcained using the ENDF/B-III parameters 
(<D,> = 6.67 eV and labeled (D) in Figure 2a) predicts consis- 
tently higher (2 to 77.) values of <SIR> than the experimental 
results. Vet the prediction of <TR> by the ENDF/B-III parameters 
is in excellent agreement with experimental <TR> as shown in 
Figure 2b. However, it is difficult to conclude that there is a 
definite parity dependence because of the possible uncertainties 
involved in the calculations (the neglect of the d-wave contri- 
bution, the theoretical ambiguity about J-dependence of strength 
functions and the experimental errors) . The experimental errors 
in <SIR> are typically +2.07. due to counting statistics and 
+ 2.0% due to estimated systematic components. 

B) Temperature Dependence of <TR> and <SIR> 

Both the ideal gas model and the Lamb effective temperature 

« 

model 11 have been applied to calculate <TR> and <SIR>, and the 
comparison is shown in Table II. For a sample temperature of 
77°K, the ideal gas model predicts a slightly larger <TR> and a 
lower <SIR> because of the larger cross-section fluctuation pre- 
dicted by this model. However, the difference between the pre- 
dictions of two models is negligible above room temperature. 

The self- shielding effect in <TR> and <SIR>, as can be seen 
in Eqs. 3 and 4, depends on: 

1) cross-section fluctuations, 

2) sample temperature, 

3) sample thickness. 

In general, the cross-section fluctuations decrease as the 
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neutron energy increases because of the decrease in peak cross 
sections, the increase in neutron width and the decrease in level 
spacing. The cross-section fluctuations also decrease as temper- 
ature increases due to the Doppler effect. Therefore, the <TR> 
and <SIR> both become smaller for lower neutron energies and for higher 
sample temperatures. These effects are also correlated with the 
thicknesses of the shielding samples. At the peak of a resonance, 
both <TR> and <SIR> increase with temperature because of the 
decrease in peak cross section at higher temperature. This 
effect is reversed at the wings of a resonance. If a shielding 
sample is thick enough to 'saturate' the transmission, in the 
vicinity of a resonance peak, then <TR> and <S1R> decrease with 
temperature; this is especially noticeable at low energies below 
2 keV. As the sample becomes thicker, the 'transmission satura- 
tion 1 effect in the vicinity of a resonance peak becomes severe, 
causing a smaller temperature effect for thicker samples than ex- 
pected. This is shown in Table III in the energy range of 200 eV 
to 300 eV for both <TR> and <SIR> of depleted uranium. As the 
neutron energy increases the 'transmission saturation' effect 
decreases so that the contribution from the resonance peak to 
<TR> and <SIR> becomes larger. This effect can be seen in 
Figure 5a where the cold-to-room <SIR> ratio decreases as neutron 
energy increases due to an increasing resonance-peak contribution 
at higher energies. 

The thermal expansion of the sample thickness also plays 
an important role in the temperature-dependence of <TR> and <SIR>. 
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Since the sample becomes effectively thinner at higher tempera- 
ture, the <TR> and <SIR> increase with temperature as can be seen 
in Table III. Here, in the energy range of 20 ‘to 30 keV, the 
observed cold-to-room and room-to-hot ratios (shown in parentheses) 
are less than unity because of thermal expansion of the samples. 

TANTALUM RESULTS 


The average resonance parameters from the analytical method 
which best fit the <TR> and <SIR> of tantalum in the energy range 
of 10 to 100 keV are listed in Table IV. In Figures 6a and 6b the 
<SIR> and < TR> based on these parameters are shown and compared 
with the experimental values for different sample thicknesses. 

The temperature dependences of <TR> and <SIR> of tantalum are 
essentially analogous to those of depleted uranium and the de- . • 

tfils* of, results have been reportedelsewhere A ^ 

The best- fit s-wave level^spacing, <D 0 >, Was 4.3^'^ eV, 
and corresponding p-wave level spacing, < D-^ > , was 2.13 ^q*|^ eV. 

The potential cross section obtained using the scattering length, 

R 1 =8.19 ferrois, ranges from 8.42 barns at 1.0 keV to 7.8 barns 
at 100 keV, both of which are slightly higher than the measured 


low energy value of 7.0 + 1.0 barns 


17 


The best- fit s-wave 


strength function (in the energy range of 10 to 100 keV) , 

S Q = 1-6^q' 2 (10 ^) > is slightly lower than the previously meas- 
ured values, which range from 1.8x10"^ (with T = 0.055 eV and 


Y 


< D o> =4.3 eV) 18 to 2.0x10'^ (with T y = 0.065 eV,<B > =4.3 eV and 

-4 19 ° 

= 0.2x10 ) . However, it appears that in the energy range 
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from 1 keV to 10 keV, a slightly higher s-wave strength function 

(S = 2.0xl0' 4 ) fits both <TR> and <SIR> better for the thinnest 
v o 

sample (Figures 6a and 6b) . 

The applicability of the analytical model depends on the 
sample thickness, the cross section fluctuations and the sample 
temperature. As can be seen in Figures 6a and 6b (and also in 
Figure 2), for the thickest sample (0.0802 atom/barn) the analy- 
tical model fails to predict <TR> and <SIR> at the energies below 
~13 keV, while the analytical model continues to work at the 
energies down to ~1 keV for the thinnest sample (0.00363 atom/barn). 
This is expected since the exponent, n-^, appearing in the ex- 
pressions of <TR> and <SIR> in Eqs. (3) and (4) becomes larger at 
lower energies and for thicker samples, and the four-term approxi- 
mation begins to break down. 

In Figure 7, the self- shielding factors calculated from <TR> 
and <SIR> , r| TR and 11 SIR > are compared in terms of the average 
effective (shielded) cross sections," 0 and * In 

general, r| TR is greater than or equal to unity, and as long as 
the potential-to-resonance interference effects of the elastic 
scattering cross section do not dominate , rig ^. R is less than or 
equal to unity that is, 




12 


Figure 7 also shows that the self- shielding effect becomes smaller 
(the value, jl-r]|, becomes smaller) as the neutron energy in- 
creases and the difference between and decreases. 


CONCLUSIONS 


This work represents an extensive investigation of the 
energy and temperature dependence of depleted uranium and tantalum 
in the unresolved resonance region. The analytical model has been 
developed for the interpretation of the data in the unresolved 
region. The results of this investigation are graphically 
summarized in Figures 2 to 7 and are listed in Tables I and IV. 

The general concluding remarks which follow from these are: 

1) The average resonance single level parameters of both 
depleted uranium and tantalum in the energy range 10 keV to 100 keV 
have been determined and compared with other measurements (Tables 

I and IV). These parameters which best fit <TR> and <SIR> agree 
with low energy parameters within the statistical counting error 
limitations for both samples, establishing the validity of the single level 
extrapolation of the resolved region parameters into the un- 
resolved resonance region. 

2) The <SIR> is more sensitive to a y and to p-wave resonance 
parameters, while the <TR> is more sensitive to a and s-wave 
resonance parameters. The combination of <SIR> and <TR> experi- 
ments provides a comprehensive determination of the key parameters 

necessary to carry out reactor calculations in the unresolved 
region. 
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3) The applicability of the analytical model can be extended 
to lower energies if the sample thickness is optimized in the 
measurements. Since the average quantities appearing in the <TR> 
and <SIR> are calculated by direct integration over the known 
statistical distribution functions, the analytical model has a 
great advantage over the statistical sampling method in calcula- 
tion time . The analytical method is 10 to 10 faster than the 
statistical sampling method. T his method extracts the aver age 
psnnarif'p paramete r s fO 1 " singln - lnwl — Br cit ■■ Wigne i foilLialiSLIlb 
-&o— that the reattTfcs eaa be directly xiaed is rmnt.nr cal enl at . 

4) The <TR> and <SIR> decrease as temperature increases in 
the energy region below ~5 keV for both the uranium and tantalum. 
However, at higher energies the <SIR> increases with temperature 
and energy, while the <TR> continues to decrease as temperature 
increases. The 'transmission saturation 1 effect at the resonance 
peak as a function of the sample thickness plays an important 
role for the temperature dependence of the <TR> and <SIR> . 

5) The ideal gas model predicts a slightly larger <TR> and 
a lower <SIR> than the predictions by the Lamb effective model. 
However, the difference in predictions between two models is 
negligible above room temperature. 
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Depleted Uranium <SIR> and <TR> Compared with the 
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resonance parameter sets labeled A, B. C and D) . 

Average Total and Capture Cross Sections of Depleted 
Uranium. 

The Sensitivity of <TR> and <SIR> to the Change of 
the Scattering Length, R ' . 

The Sensitivity of <SIR> and <TR> to the Change of 
the Strength Functions. 

The Sensitivity of <SIR> and <TR> to the Change of 
the Level Spacings. 

The Cold-to-Room Ratio of <TR> and <SIR> for Depleted 
Uranium of 0.03155 and 0.06206 atom/barn. 

Tantalum <SIR> and <TR> Compared with the Analytical 
Calculations . 

Effective Average Cross Sections of Tantalum Calculated 
from <SIR> and <TR>. 
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The Best-Fit Average Resonance Parameters of Tantalum. 



Table I. Best-Fit p-wave Strength Functions for 


238 U for Different Scattering Lengths 



Best Fit Parameters* 

ENDF/B-III 

(D) 

(E) 

(A) ] (B) 

(C) 

SjUO' 4 ) 

1 SR 

1 ” 58 -0.1 

1 94 d 

L - y4 -0.2 

2 - 40 -0.2 

1.75 

1.27 

<D t > (eV) 

11 3 “ 2 -° 

li. a +2 „o 

11 3 +3 - 5 

11 +3.0 

u • 3 il:° 0 

6.67 

6.90 

f Yl (meV) 

47 ‘ 5 +8.4 

43 8 “ 13 ’ 6 
43 - 8 +11.6 

37.0 

1 

28.0 

R ' (fin) 

9.30 

j 9.20 

9.0 

9.20 

9.60 

R (fm) 

8.74 

8.74 

8.4 

8.4 

8.74 

R 90 

-0.065 

-0.053 

-0.071 

-0.095 

-0.098 

s o ( 10 - 4 ) 

.. ...i ■■ ■ 

1.0 

1.05 

1.0 

<D q > (eV) 

20.7 

20.0 

20.7 

f (meV) 

yo 

23.0 

23.5 

23.0 


*The uncertainties given for S l5 <V l >, and r 
difference between the experimental values, 



N 

I \ 

N L, 
i=l 




)/* 



2,0 


correspond to the values which yield the average 
r? x P , and theoretical values, , as follows: 


where 6R? Xp denotes 


the statistical counting errors 


of the results. 



Table II. Depleted-Uranium <TR> and <SIR> Calculated from Ideal 
Gas Model and from the Effective Temperature The or/ 


Energy 

Range - 30.6 - 27.7 V.eV 


S amp 1 e 

Temperature = 77"’K 


Sample 

Thickness = 0.03155 at/'o 



"TR> 

; - 

Experimental 

0.658 

r 

Values 

-0 .00 7 

y • ^ - 

Effective Temp. 
Theory 

C . 6f 9 

<r ~ • 

Ideal Gas 
Model 

0.661 

- - V 



Table III. 


Cold-to-Room (and Room- to-Hot) Ratios of Depleted -Uranium <TR> and <SIR 
in Typical Energy Regions 


Energy 

Regions 

(keV) 

Sample 

Thickness 

(at/barn) 

"/C 

Cold-to-Room Ratio 

Room-to-Hot Ratio * 

<TR> 

<SIR> 

<TR> 

<SIR> 


0.03155 

1.019(4) # 
(1.016) + 

1.215(9) 

(1.212) 

1.026(4) 

(1.006) 

1.293(8) 

(1.268) 

0.2-0. 3 

0.06206 

1.008(3) 

(1.003) 

1.076(5) 

(1.071) 

1.030(4) 

(0.989) 

1.350(12) 

(1.294) 


0.03155 

1.017(3) 

(1,015) 

1.077(5) 

(1.075) 

1.035(3) 

(1.017) 

0.986(4) 

(0.969) 

1. 0-2.0 

0.06206 

1,021(3) 

(1.017) 

1.078(3) 

(1.073) 

1.048(2) 

(1.014) 

1.094(4) 

(1.055) 

20.-30, 

0.03155 

1,005(3) 

(1.003) 

0.995(5) 

(0.993) 

1.007(3) 

(0.997) 

0.977(5) 

(0.962) 

0.06206 

0.996(6) ‘ , 

(0.993) 

— — 

0.985(4) 

(0.982) 

1.006(3) 

(0.975) 

0.993(6) 
. (0.962) 


" Sample temperatures are 77 , 295 and 973 K. 


# 


+ 


The number, 1.019(4), 


:TR> 


:old 


/ <TR: 


room 


and this value is for 
The values inside the 


stands for the ratio, 

= 1.019 + 0.004 (statistical error only) 
the ratio after the samples are corrected for 
parentheses are the ratios before the sample 


thermal expansion, 
thickness corrections. 



Table IV. Best-Fit Average Resonance Parameters of Tantalum 
in the Energy Range from 10 to 100 keV 



The error limitations given for all the parameters correspond 
to the values which yield the average difference between the 
experimental values, R? Xp , and theoretical values, R? , as follows 
N 1 ' 

AK = ^ r _ (R l Xp - R^ al )/«Ri Xp £ 2.5 

i=l 

where SR? xp denotes the statistical counting errors of the results. 

It is assumed that the channel radius R is the same as R ' , 
that is . R 00 = 0.0. 


